Magnetoresistive junctions based on epitaxial graphene and hexagonal boron nitride 
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We propose monolayer epitaxial graphene and hexagonal boron nitride (ft-BN) as ultimate thick- 
ness covalent spacers for magnetoresistive junctions. Using a first-principles approach, we investigate 
the structural, magnetic and spin transport properties of such junctions based on structurally well 
defined interfaces with (111) fee or (0001) hep ferromagnetic transition metals. We find low resis- 
tance area products, strong exchange couplings across the interface, and magnetoresistance ratios 
exceeding 100% for certain chemical compositions. These properties can be fine tuned, making the 
proposed junctions attractive for nanoscale spintronics applications. 

PACS numbers: 72.25.-b, 73.43.Qt, 75.47.-m, 81.05. Uw 
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I. INTRODUCTION 

Graphene, a recently discovered two-dimensional form 
of carbon, has attracted unrivaled attention due to its 
unique physical properties and potential applications in 
electronics i^i^ This nanomaterial is particularly promis- 
ing for the field of spintronics, which exploits both the 
spin and the charge of electrons .^^'^'^'^'^1^ One funda- 
mental spintronic effect is the magnetoresistance, the 
change in electric resistance as a function of the rel- 
ative orientation, either parallel or antiparallel, of the 
magnetization of two ferromagnetic layers separated by 
a nonmagnetic spacer layer fiSi Achieving high magnetore- 
sistance ratios while keeping reasonably low electric re- 
sistance is crucial for many technological applications.^- 
However, reaching this goal is currently hindered by 
material-specific restrictions such as the inability of pro- 
ducing well-ordered ferromagnet/spacer interfacesJ^ii^ 

Semimetallic graphene and its insulating counter- 
part, isostructural hexagonal boron nitride (/i-BN), 
are promising spacers as epitaxial monolayers of these 
materials can be grown by means of chemical va- 
por deposition (CVD) on a broad variety of metallic 
substrates] The quality of such epitaxial 
monolayers is very high and the covalent bonding net- 
work of both graphene and /i-BN is perfectly preserved 
upon the bonding to the substrate. Moreover, the growth 
of graphene and /i-BN on fcc(lll) and hcp(OOOl) sur- 
faces of ferromagnetic Co and Ni results in commensu- 
rate epitaxial layers due to the closely matching lattice 
constantsjii This has led to a theoretical prediction of 
perfect spin filtering and, thus, to extremely high magne- 
toresistance ratios in such junctions based on multilayer 
graphene (>4 layers) and graphite. However, the CVD 
growth on crystalline surfaces is self-inhibiting, that is 
only one epitaxial layer can be grown. The deposition of 
ferromagnetic nanoparticles on top of epitaxial /i-BN has 
also been demonstratedi ^^'^^ These interfaces further of- 
fer the opportunity of fine tuning their properties through 
the intercalation of other metals, such as FcfSi Cr^ and 
Au^ 

In this work, we suggest the use of monolayer graphene 



and /i-BN as covalently bonded spacer layers of minimal 
thickness in magnetoresistive junctions. Through first 
principles calculations we study the structural, magnetic 
and spin transport properties of such junctions based 
on first-row ferromagnetic transition metals: natural hep 
and fee Co, fee Ni, as well as intercalated fee Fe. We show 
that the proposed magnetoresistive junctions realize low 
electric resistances, strong interlayer exchange couplings, 
and magnetoresistance ratios exceeding 100% for certain 
chemical compositions. 

This paper is organized as follows. In Sec. |TT] we 
describe our computational methodology, including the 
first-principles approach to electronic transport. In Sec. 
mil we report the atomic structure and electronic proper- 
ties of the considered magnetoresistive junctions. Partic- 
ular attention is devoted to the interlayer exchange cou- 
plings. The results of electronic transport calculations 
are discussed in Sec. IIVI Section |V] concludes our work. 



II. COMPUTATIONAL METHODS 

The electronic and atomic structure calculations were 
performed using the pwscf plane-wave pseudopoten- 
tial code of the quantum-espresso distribution.^^ To 
achieve a good description of atomic structures, in- 
terlayer exchange couplings and spin transport prop- 
erties, we chose the Perdew-Burke-Ernzerhof exchange- 
correlation density functional. Ultrasoft pseudopoten- 
tials were used to describe core-valence interactions 
The valence wave functions and the electron density were 
described by plane-wave basis sets with kinetic energy 
cutoffs of 25 Ry and 250 Ry, respectively.— The atomic 
structure of the magnetoresistive junctions considered in 
our work is illustrated in Fig. [T]^a). Our investigation is 
restricted to only symmetric junctions, i.e. with the same 
metal on both sides of the spacer layer. Each ferromag- 
netic layer consisted of six atomic planes. The solutions 
for parallel and antiparallel relative spin orientations of 
these two layers were obtained by specifying appropriate 
initial orientations of the magnetic moments. The lateral 
unit cell of the studied interfaces is shown in Fig. [^b). 
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FIG. 1: (Color online) (a) Representation of the atomic struc- 
ture of magnetoresistive junctions based on epitaxial mono- 
layer graphene and /i-BN. (b) Top-view of graphene on (111) 
surface of fee Co or Ni. The two-dimensional unit cell is indi- 
cated by dotted lines and the principal atomic positions are 
labeled, (c)-(e) Side-views along the longest unit cell diago- 
nal of the lowest-energy interfaces formed by Co (hep and fee) 
and Ni (fee) in combination with either monolayer graphene 
or /i-BN. 



We considered bound configurations and determined the 
lowest-energy structures through the relaxation of atomic 
positions. For these configurations, we performed quan- 
tum transport calculations in the current perpendicular 
to plane configuration using the pwcond code'^° of the 
same package. The scattering region included the spacer 
monolayer and three adjacent monolayers of metal on 
both sides. We use the optimistic definition of the mag- 
netoresistance ratio: 
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The spin-resolved quantum conductances G^ for paral- 
lel (cr =tT,J,J, for majority and minority spins, respec- 
tively) and antiparallel {a =ti) configurations were cal- 
culated by integrating the corresponding k|| -dependent 
transmission probabilities T^^ evaluated on a uniform 
grid of 64x64 k-points in the two-dimensional Brillouin 
zone. 



III. ATOMIC AND ELECTRONIC STRUCTURE 

To determine the lowest energy structures of the junc- 
tions, we carried out structural relaxations for all pos- 
sible stacking orders of the atomic planes in the vicin- 
ity of the spacer layer. The corresponding structures 
are shown in Fig. [T] for Co and Ni based junctions and 
summarized in Table U for all investigated chemical com- 
positions. We find that both graphene (GR) and h-BN 
bound to the transition metals (TM) display short metal- 
carbon and metal-nitrogen distances (2.19-2.45 A) com- 
parable to the sum of the corresponding covalent bond 
radii. The thickness of the spacer layer is thus compara- 
ble to that of a single atomic plane of the ferromagnetic 
metal. For the Ni|GR|Ni(fcc) junction we find a Ni-C 
distance of 2.19 A, which is close to 2.18 A calculated for 
the graphene chemisorbed on the Ni(lll) surface (i.e. 
Ni|GR system). The latter value is in good agreement 
with the experimental value of 2.16±0.07 A.^"* 

We now turn to the electronic structure of these mag- 
netoresistive junctions. Figure [2] shows the spin- resolved 
projected density of states (PDOS) onto the light atoms 
(B, C and N) of the fee Co junctions. One can see that 
the characteristic "Dirac cone" density of states of the 
free-standing graphene is not preserved upon the forma- 
tion of the Co|GR|Co(fcc) interface [Fig. [^a)] . This is 
consistent with the theoretically predicted^" and exper- 
imentally observedii strong hybridization between the 
electronic states of graphene and of the TM surface. Sim- 
ilarly, both B and N centered states fill the band gap of 
the insulating /i-BN [Fig. [2Ib)]. In both cases, we very 
similarly find significant contributions of the epitaxial 
layer states to the density of states at the Fermi level. In 
the parallel (antiparallel) configuration of the graphene 
based junction, the induced magnetic moments on the 
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FIG. 2: (Color online) Spin-resolved projected density of 
states (PDOS) onto the atoms of the spacer layer, either (a) 
graphene or (b) /i-BN, in the fee Co junctions in their parallel 
(P) and antiparallel (AP) configurations. The majority and 
minority spin labels refer to the parallel configuration; in the 
antiparallel arrangement the spin channels are equivalent. 
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FIG. 3: (Color online) k|| -Resolved conductance per unit cell 
(in units of /h) through (a) bulk hep Co along the (0001) 
direction, (b) Co|GR|Co(hcp) and (c) Co|BN|Co(hcp) junc- 
tions, and (d) a vacuum layer of equivalent thickness. The 
columns correspond to majority and minority spin channels 
of the parallel configuration, and to one of the equivalent spin 
channels of the antiparallel configuration, respectively. Labels 
indicate the total conductances per unit cell area. 



carbon atoms in the unit cell are — 0.005/iB (0.081/iB 
and — 0.081/xb). In the parallel configuration of the h- 
BN junction, the induced magnetic moments of N and 
B atoms are 0.029/iB and — 0.065/iBi respectively. In the 
antiparallel configuration both vanish by symmetry. 

The interlayer exchange coupling, the difference Ai? = 
Ep~Eap between the energies of parallel and antiparallel 
configurations, is a manifestation of the superexchange 
mechanism. It achieves rather high values [cf. Table [I] 
due to the ultimate thickness of the spacer layer. For 
Fe and Co, the antiparallel configuration is energetically 
favored. On the contrary, the parallel configuration is 
preferred for Ni. This intriguing crossover provides op- 
portunity for fine tuning the interlayer exchange by vary- 
ing the chemical composition of the ferromagnetic layers. 



IV. ELECTRONIC TRANSPORT 

A. Role of spacer material 

To understand the calculated quantum conductances 
and the resulting magnetoresistance ratios [cf. Table |T] , 



we analyzed the k|| -resolved transmission probabilities. 
First, we studied the effect of the spacer layer in hep Co 
junctions which have the same lowest energy structure 
for both graphene and /i-BN [Fig. [3]. We found that 
both systems show strikingly similar k|| -resolved trans- 
mission probability maps [compare Figs. [3Jb) and[3]Jc)] 
and consequently quantum conductances 



T, 



The tIJ^ and 



1^1 1 maps reveal major features of the projected hep Co 
Fermi surfaces for the free-electron-like majority spin and 
mostly d-symmetry minority spin electrons [Fig. OJa)], 
which are relevant to the quantum conductances of bulk 
metals. ■^^'•^^ The total transmission probabilities of the 
junctions in the parallel configurations constitute ~40% 
and ~20% of the Sharvin conductances'^'* of bulk hep Co 
along the (0001) direction. The quantum conductances 
in the antiparallel configuration are mostly determined 

Their values are conse- 



by the overlap of rj|j 



and 

0.140 e'^/h per unit cell for both 



quently lower {G-\i 
spacer materials). The resulting magnetoresistance ra- 
tios are 86% and 66% for graphene and /i-BN junctions. 
Thus, in the regime of ultimate thickness the transport 
properties are largely independent of the electronic struc- 
ture differences of the two spacer materials. The role 
of a single layer of covalent spacer material consists in 
fixing a certain stacking order at the interface and in 
providing a medium for the abrupt change of magnetiza- 
tion in the antiparallel configuration. Due to the metal- 
lic nature of the spacer layers [cf. Fig. [2] such junctions 
possess low resistance area products (<3xl0^^^ flui^) 
which makes them suitable for nanoscale spintronics ap- 
plications such as the magnetic random access memo- 
ries and spin transfer nano-oscillators. We classify the 
present systems as giant magnetoresistance (OMR) junc- 
tions. This contrasts to the spin transport through a 
vacuum gap of the same thickness which shows de- 
caying with |k||| [Fig. El^d)], a characteristic feature of 
tunneling. ■^^ ^j^g magnetoresistance ratio is about twice 
smaller (38%) in the case of tunneling through a vacuum 
gap. This allows us to conclude that in the limit of ulti- 
mate thickness the GMR effect is more efficient than the 
tunneling magnetoresistance. 



B. Role of ferromagnetic layers 

Next, we studied the dependence of transport prop- 
erties on the ferromagnetic metal by considering fee Fe, 
Co and Ni junctions in combination with /i-BN. For all 
three metals, the majority spin transmission in the paral- 
lel configuration, rj|^, undergoes little change along the 
Fe-Co-Ni series [Fig. 2]. This behavior stems from the 
similarity of the corresponding majority spin Fermi sur- 
faces of the bulk metals, which are formed by partially 
filled s bands. However, much larger differences are found 
for Tj^ij involving the minority electrons. These reffect the 
drastically different Fermi surfaces resulting from the in- 
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TABLE I: Lowest energy structures, interlayer exchange couplings {AE — Ep — Eap), spin-resolved quantum conductances 
(G||, Gil ^-iid Gil) s-nd magnetoresistance ratios (MR) for the discussed graphene (GR) and /i-BN junctions. The notation 
for the stacking order is identical to the one in Fig. [1] The values of quantum conductances and interlayer exchange couplings 
are given per unit cell area. 



junction 


stacking 


AE 


G|l 


Gi 


Gil 


MR 




order 


(meV) 






(e^h) 


(%) 


Fe|GR|Fe(fcc) 


cba\bac 


79 


0.334 


0.440 


0.240 


61 


Fe|BN|Fe(fcc) 


cba\abc 


63 


0.256 


0.297 


0.111 


149 


Co|GR|Co(fcc) 


bca\bca 


91 


0.317 


0.427 


0.232 


60 


CoBN|Co(fcc) 


bca\acb 


46 


0.263 


0.268 


0.210 


26 


Ni|GR|Ni(fcc) 


bca\bca 


-18 


0.352 


0.587 


0.402 


17 


NiBN|Ni(fcc) 


bca\acb 


-3 


0.207 


0.722 


0.299 


55 


Co|GR|Co(hcp) 


aca\aca 


29 


0.241 


0.278 


0.140 


86 


CoBN|Co(hcp) 


aca\aca 


44 


0.222 


0.241 


0.140 


66 



terplay between s and d states. The increase of Gn along 
the series can be attributed to the decrease of hybridiza- 
tion between s and d electrons upon the increase of d 
band fiUingj^^ in general, the free-electron-like s states 
show higher transmission probabilities. The Gn values 
are again determined by the overlap of T^J^ and T^j'^ and 
tend to increase along the series. For the Fe|BN|Fe(fcc) 
junction we find a magnetoresistance ratio of 150%, the 
largest value among the compositions studied. 

Further search of magnetoresistive junctions with im- 
proved characteristics may consist in exploring asymmet- 
ric junctions and the intercalation of some other chemical 
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FIG. 4: (Color online) k||-Resolved conductance per unit cell 
(in units of e^/h) through fee (a) Fe, (b) Co and (c) Ni junc- 
tions based on monolayer /i-BN. The columns correspond to 
majority and minority spin channels of the parallel config- 
uration, and to one of the equivalent spin channels of the 
antiparallel configuration, respectively. Labels indicate the 
corresponding total conductances per unit cell area. 



elements at the interfaces. We here demonstrate the sec- 
ond possibility. It has been suggested that the incorpora- 
tion of submonolayer quantities of Cu at the TM|GR„ in- 
terface would reduce undesired hybridization between the 
states of graphene and of the metal surface at the price of 
substantially decreasing the magnetoresistance ratiosi^ 
However, we find that the decoupling of the spacer layer 
from the metal surface does not necessarily imply the 
loss of magnetoresistance. This can be achieved by in- 
tercalating the metals from the middle of the transition 
metals series, e.g. Mn, which show reduced binding to 
carbon tt systemsj^ Indeed, in the intercalated CoMn(l 
ML)|GR|Mn(l ML)Co(hcp) junction the Mn-C distance 
increases to 2.95 A and the interlayer exchange coupling 
decreases to 10 meV (to be compared with 29 meV for 
Co|GR|Co(hcp), cf. Table IJ). Concurrently, the magne- 
toresistance ratio raises from 86% to 127%. The Mn layer 
is strongly spin polarized and antiferromagnetically cou- 
pled to hep Co. 



CONCLUSIONS 



In conclusion, we propose epitaxially grown monolayer 
graphene and /i-BN as ultimate thickness covalent spac- 
ers in transition metal based magnetoresistive junctions. 
Such junctions display well-ordered interfaces and can 
be produced through existing manufacturing processes. 
Their physical properties can be fine tuned in a broad 
range by varying the chemical composition. These sys- 
tems show low resistance area products and typical GMR 
behavior with magnetoresistance ratios exceeding 100% 
for certain compositions. Both ferromagnetic and anti- 
ferromagnetic interlayer exchange couplings are found. 
These properties make the proposed junctions attractive 
for spintronics applications such as the magnetic random 
access memories and spin transfer nano-oscillators. 
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